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We introduce a formalism for the calculation of the time of arrival t at a space point for particles 
traveling through interacting media. We develop a general formulation that employs quantum 
canonical transformations from the free to the interacting cases to construct t in the context of the 
Positive Operator Valued Measures. We then compute the probability distribution in the times of 
arrival at a point for particles that have undergone reflection, transmission or tunneling off finite 
potential barriers. For narrow Gaussian initial wave packets we obtain multimodal time distributions 
of the reflected packets and a combination of the Hartman effect with unexpected retardation in 
tunneling. We also employ explicitly our formalism to deal with arrivals in the interaction region 
for the step and linear potentials. 

1 PACS numbers: 03.65.Bz, 03.65.Ca, 03.65.Nk 
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I. INTRODUCTION 

In this paper we work out a theoretical framework to compute the time in which a particle that moves through an 
interacting medium arrives at a given point. In the construction of this framework we will have to deal with problems 
of very different kind that we introduce now: 
^| ' First, there is the nature of time in quantum mechanics. It appears as the external evolution parameter in the 
Schrodinger and Heisenberg equations, common to both, systems and observers alike. However, time arises in many 
■ instances (transitions, decays, arrivals, etc.) as a property of the physical systems. The attempts to promote time to 
the category of observable run early into the obstruction detected by Pauli Q: A self-adjoint time operator implies 
r"| ■ an unbounded energy spectrum. This was soon related to the uncertainty relation for time and energy, whose status 
Q H | and physical meaning produced some controversies [||-|| , and is still subject of elucidation today (see for instance || 
and J/J). The question remains unsettled for closed quantum systems, specially in the case of quantum gravity, whose 
' formulation is pervaded by the so called problem of time || . 

Second, the definition of the time-of-arrival (toa), which is probably the simplest candidate time to become a 
property of the (arriving) physical system, rather than a mere external parameter. Due to its conceptual simplicity, 
it has been used in many cases to illustrate different problems related to the role of time in quantum theory. Allcock 
analyzed |^] extensively the difficulties met by the toa concluding that they were insurmountable. The present 
situation is ambiguous. On the one side, there are theoretical analysis fl(]Jll| of the toa suggesting that it can not 
be precisely defined and measured in quantum mechanics. This contradicts the possibility of devising high efficiency 
absorbers |l2|, that could be used as almost ideal detectors for toa |L3]]. On the other hand, there are explicit 
constructions of a self-adjoint (albeit in a pre-Hilbert space) toa operator for the non-relativistic free particle in one 
space dimension |TJ,|l^|, and for the relativistic free particle in 3-D p6[ , both avoiding the Pauli problem. There 
is also an alternative formulation Jl7j as a Positive Operator Valued Measure (POVM). Finally, the toa has been 
measured in high precision experiments Jl8| , ^9[ on the arrival of two entangled photons produced by parametric down- 
conversion, one of which has undergone tunneling through a photonic band gap (PBG). The experimental results that 
show superluminal tunneling, neatly identify the Hartman effect ^lj and the Wigner time delay |^] (or phase time) 
as the physically relevant mechanisms for the tunneling time and toa respectively. Whether these results apply only 
to photons and are due to the specific properties of the PBG used, or can be extended to other particles and barriers, 
can not be decided in the lack of a satisfactory theory of the toa at a space point through interacting media. 

The third question is thus the tunneling time, for which there are three main proposals. Wigner introduced the 
phase time in his analysis (22) of the relationship between retardation, interaction range, and scattering phase shifts. 
Buttiker and Landauer introduced the traversal time |23| in their study of tunneling through a time-dependent barrier. 
Soon after, Buttiker used the Larmor precession as a clock [Q, identifying the dwell [2q] , traversal, and reflection 
times as three characteristic times describing the interaction of particles with a barrier. Recent reviews that include 
these and other approaches, discussing toa and tunneling times from a modern, unified perspective, can be found 
in |2q] and [E7| . The light shed on these questions by the two photon experiments is revised in [Eq] and |23| . 
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The main progress, quoted before, towards the formulation of a quantum toa operator has been its explicit con- 
struction for a particle moving freely in one space dimension. In this paper we face the problem of extending this 
formalism to the case of the presence of an interaction potential affecting a region of the (1-D) space. We extend 
here our previous unpublished results f3(i|| , taking now proper care of the dependence on the arrival position x that 
we consider placed in front of the interaction region, within it, and behind it (if spatially finite). 

The plan of the work is as follows: In Section II we construct a toa formalism of general validity. The starting point 
is the case of the free particle. There, a suitable canonical transformation, the quantum version of the Jacobi-Lie 
transformation of classical mechanics, gives the toa in interacting media, (even at points where V(x) ^ 0). In Section 
III we consider an initial state consisting of a narrow Gaussian wave packet prepared at the left of the interaction 
region and moving towards it. A quasi-classical study of the toa at a point x in the interaction region is first carried 
out. Then we turn to the full quantum mechanical treatment. We first analyze the arrival in the presence of a step 
potential. Section IV is devoted to the study of the toa at points behind square barriers. We detect, in diferent 
instances, stauration of and departures from the Hartman effect. The case of x at the left of the interacting region 
is characterized by the (possibly interfering) contributions coming from the incident and the reflected wave packets. 
This situation is treated in Section V, where we deal separately with the case of total reflection (very high barriers) 
which has an analog in classical mechanics, and with the case of partial reflection, a pure quantum phenomenon with 
very rich structure in the time domain. Finally, we summarize our results in Section VI. 

In Appendix A we show how the toa can be treated as a derived quantity in the phase space of classical Hamiltonian 
systems in the case where these are integrable. A short review of the construction and properties of the quantum toa 
operator for free particles is presented in Appendix B for completness. 



To measure the time of arrival of a free particle at a point x one would: a) place a detector at x, b) prepare the 
initial state |?/>) of the particle at t = 0, and then, c) record with a clock the time t when the detector clicks. The 
value of t gives the toa of the state \ip) at x. Repeating this procedure with identically prepared initial states, one 
would get the probability distribution in times of arrival at x. Of course, the results would depend on the initial state 
chosen, which stores all the information regarding the initial distribution in positions and momenta of the particle. 

We want to determine the effect on these times of a position dependent interaction between the particle and the 
medium, that we describe by a potential energy V(q). For instance, to disclose the effect of climbing (or tunneling 
through) a potential barrier, one would simply put the barrier in between the detector and the initial state, and then 
record the new times of arrival. With an initial state identical to that prepared for the free case, any difference in 
the probability distributions should be an effect of the barrier. Several questions can be investigated by changing the 
properties of the barrier: its height or width if it is rectangular, even its very form. This has been explicitly done in 
the two photon experiments at Berkeley, by putting alternatively a mirror and an ordinary glass in the path of one of 
the photons. Of great interest is the dependence on x where V(x) > 0, i.e. with the detector within the range of the 
interaction, and also the time of arrival for E > V, in which classically there is no reflection. Some of these questions 
are studied in last sections of this paper. 

In classical mechanics particles move along the trajectories H(q,p) =const. as t increases. This allows to work out 
t x , the time of arrival at the point q(t) = x, by identifying the point (q,p) of phase space where the particle is at (say) 
t = 0, and then by following the trajectory that passes by it, up to the arrival at x. The mathematical translation of 
this procedure is given by the equation of time: 



that is discussed at length in many textbooks, and whose existence conditions and characterization as a function of 
the phase space variables are outlined in the Appendix A. We simply note here that t x (q,p) is canonically conjugate 
to the hamiltonian {t x (q,p), H(q,p)} = —1. 

This equation is a troublesome starting point for quantization. First, it involves a (path) integral of operators and 
should be treated accordingly. Second, it only applies to values of x that are classically within the reach of (<?,p), 
while in quantum mechanics all values of x are attainable. Classically, the particle propagates without reflection up 
to the turning point qo (V(qo) = E), where it is completely reflected. There is no further penetration beyond this 
point. The situation is different in quantum mechanics: there may be tunneling beyond qo and partial reflection before 
reaching it. These phenomena cannot be accounted for by Eq. (|l|), that gives complex numbers for these cases. Now, 
note that both, tunneling and partial reflection, are absent from the motion of free particles, whose time of arrival 
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has been successfully quantized as said in the introduction. In addition, all the positions are within the reach of the 
free particle. Summarizing, everything points to the free time of arrival as a main clue to solve the problem. 

In this work we desist from attempting the straightforward quantization of the classical expression (|]) . Instead, we 
will construct the solution to the interacting case taking as starting point the well known results that apply to the 
free case. The aim is to produce the quantum version of the Lie transformation from the actual flow in phase space 
to the canonically equivalent parallel flow of constant velocity translations. In other words, we shall use the quantum 
version of the canonical transformation to action-angle variables. The Lie procedure - that we sketch for completncss 
in Appendix A - has a property that will be the crux of the matter in our construction. Namely, it permits to 
define time as a derived variable in phase space in terms of the free action-angle variables as well as, alternative and 
equivalently, in terms of the original positions and momenta. Obviously, both definitions give the same result as we 
show explicitly in Eq. ( |6l| ) . Our use of the Lie procedure in the quantum case can be described as the combination 
of steps a,b and c below: 

a. The quantization of the time of arrival to of the free particle. This is an old problem in quantum mechanics, 
whose solution in terms of a Positive Operator Valued Measure we describe in Appendix B. 

b. The construction of the quantum canonical transformation U that connects the free particle dynamics with 
Hamiltonian Hq to the case of interest with Hamiltonian Hq + V(q). U is given by the Moller wave operator as 
we show in Sections 2.1 and 2.2. 

c. The application of the canonical transformation U to to to get the time of arrival t in the presence of the 
interaction potential V(q), that is t = U to U*. This is what we do in Section 2.2, where we also address the 
interpretation of the resulting formalism. 



A. Implicit quantum canonical transformations 

Classical canonical transformations q — q(q,p), p = p(q,p) in phase space can be defined implicitly by the use of 
auxiliary functions F, G, F, G in the following way: 

F(q,p)=F(q,p) 

G(q,p) = G(q,p). (2) 

It is easy to work out the following relation among Poisson brackets: 

{F,G} qp {q,p} qp = {F,G} qp (3) 

In these conditions, the transformation is canonical (i.e. {q, p} q p = 1) if and only if 

{F,G} qP = {F,G} qp . (4) 

This relation has the additional property of fixing one of the four functions F, G, F, G, once the other three are given. 
We can choose F and G as the free particle Hamiltonian and time of arrival respectively. Then, if F is the complete 
Hamiltonian H, G will be the corresponding toa t x given by (|l|) along the classical trajectories. 

Canonical transformations were introduced by Dirac in quantum mechanics |3l|] by the use of unitary transforma- 
tions U {UU^ = U^U = 1 )• If the operators q, p are canonically transformed from q, p , then there is a unitary 
transformation U such that 

q = U^qU 

p = U^pU. (5) 

Then one can define implicitly quantum canonical transformations, like the classical ones. This possibility has been 
thoroughly analyzed and developed. The main results of the method are collected in p2] |, where one can also find 
references to other relevant literature. The transformation U is given by 



F(q,p) = E/tF( q , p )[/ = F(q,p) 
G(q,p) = C/ t G(q,p)C/ = G(q, P ; 



(6) 

where the last equality in each row is the definition of the barred operators in terms of the unbarred ones, while the 
first equality comes from the straight application of (|^) to the l.h.s. Being U a unitary transformation, the spectra of 
the canonically transformed operators have to coincide, that is: 
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cr(q) = o-(q) = ft, cr(p) = cr(p) = ft 

<7(F)= ( r(F), ct (G)=<t(G) (7) 

where the second row stands because F and F, G and G are also unitarily related operators. 

The above relations permit to build the operator G once F, G and F are given. We assume that F and F are 
self-adjoint operators, with the eigenstates corresponding to the same eigenvalue A/ given by: 

F|/)=A / |/), F|/)=A/|/) (8) 

They form orthogonal and complete bases satisfying 

(M/V) = 8 ss ,5(\ f - X' f ), £ s J a{x) dX f \f s )(fs\= I (9) 

(M/V) - S ss ,5(X f - X' f ), J2 S L(\) dX f \fs){fs\= I (10) 

where we allow for some degeneracy (that has to be the same for both F and F) labeled by s. We have also assumed 
that A is continuous, while s is a discrete index. These assumptions could be changed straightforwardly if it were 
necessary. Now, an operator U satisfying the first row of Eq.(^) can be given simply as: 

C/ = E/ dX f\f~ s )(f s \ (u) 

It is straightforward to verify that it is unitary. We can now proceed to the sought for result: the definition of G in 
terms of G using U, that is G = UGU^. The full fledged expression is 

G(q,p) = W dX f dX f ,\f S )(f S \G(q,p)\fs>)(fi S '\ (12) 
that constitutes our main result in the quantum canonical formalism. 



B. Definition of the time of arrival 

We will now apply the above to the case where F is the free Hamiltonian Hq, F the complete Hamiltonian H and 
G the time of arrival of the free particle Eq. (|62|) . Then, we have Hq = U^HU and ITo(:r) = U^H(x)U . In Appendix 
B we have summarized the toa formalism for the free particle, given by the positive operator valued measure Pq of 
Eq.(^37j). Accordingly, the POVM P of the interacting case will be given by (c/(Q)) 

P(n{x);t u t 2 ) = UP {n (x);t u t 2 )Ul (13) 

Finally, the time of arrival operator in the presence of interactions (the G of our problem) is given by 

t(H,U(x)) = Uto(HoMx))tf. (14) 

Three comments are in order here: 
• Fixing U by the relation between both hamiltonians leads to two different solutions: 

/>oo 

U {±) = V / dE\Es(±)){EsO\ = n (±) (15) 
s Jo 

which are the Moller operators connecting the Hilbert space Tti n and TC ou t of free particle states to the Hilbert 
space Tt of the bound and scattering states. These operators are only isometric in the presence of bound states, 
because the correspondence between states in TC and free states can not be one to one. In this paper we will 
consider only well behaved potentials (V(q) > Vg), that vanish at the spatial infinity, for which the Moller 
operators are unitary because there is one free state for each scattering state. In this case, the intertwining 
relations Htt± = £1±Hq can be put in the usual form H = fl±Ho£l±. In addition, we shall adhere to the standard 
conventions, choosing (with E — \im e ^ +(E + ie)) in (|Ilf ) that gives signal propagation forward in time. 
The results that would be obtained with ^(-) would correspond to the time reversal of the actual situation. If 
t is the time reversal operator P(_)(ll(a;); t\, t 2 ) — r P(+) (n(a;); ~t 2 , — 1\) . For notational simplicity, we will 
omit this label (+) wherever possible. 
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• The reduction of the problem to a sort of free particle problem by means of a canonical transformation as done 
in ( fbt ) should not be a surprise. On the contrary, this is the quantum counterpart of the classical situation 
where the trajectories of completely integrable phase space flows can be straightened out to those of a free 
particle by means of a canonical transformation. To the classical Lie transformation of Appendix A that carries 
out this stretching corresponds the quantum transformation described above and in the previous Section 2.1. 
Concretely, Equation ( |6l| ) is the classic analog to (|l4|). 

• x is the actual detector position in the interacting case. Therefore, the arguments of t in (|l4|) have to be 
LT(a;) = \x)(x\ and H. This gives for the argument of to an object IIo(a;) = £tfTl(x)£l which is not a position 
projector. Instead, it collects all the states of the free particle that add up to produce the position eigenstate \x) 
of the interacting case by the canonical transformation. Much of the difference between the classic and quantum 
cases is hidden here, in particular the quantum capability to undergo classically forbidden jumps in phase space. 

Summarizing, in the interacting case we have a toa operator given by 

/+oo 
dtt\txs)(txs\, (16) 

where 

\txs) = {—) 1 ^e iHt Tl a \x) (17) 
m 

Above we have introduced the projector U s = J dE \Es(+))(Es(+)\, which is obtained from the II s o of Eq. ( |66| ) by 
the canonical transformation (|l5|). We now have the tools necessary for a physical interpretation in terms of a POVM: 
Given an arbitrary state ip at t = 0, its time of arrival at a position x has to be, according to (^6|), 

1 r +oc 

Mt*M = ^ Ey ^ dtt\(tx S m\ (is) 

with the standard interpretation of J2 S l(^ xs IV')| 2 as the (yet unnormalized) probability density that the state 
arrives at x in the time t. The probability of arriving at x at any time is then P(x) = J dt ^ s |(tes|V')| 2 , giving a 
normalized probability density in times of arrival 

normalization that has been used in (|l8|). Note that in the cases where P(x) vanishes this conditional probability is 
devoid of meaning: If there are no arrivals at all, there are no arrivals in any finite (or infinitesimal) interval of time. 

The above equations ( jl8] , |l9| ) can be given a form that is very useful for computation, while throws some light on 
the physical meaning of the different quantities involved. By using explicitly ([l7|), one gets 

P(x)=J2i [ dE(—)^{x\Es(+)){Es(+mr x 
— ' J m 

{ J dE 1 {^-f'Sx^' s{+)){E' s(+)\^)} J dte-« E - E '» 

= 2nJ2 I dE( — )^\(x\E S (+))(E S (+M)f (20) 
J m 



Using a similar procedure, one gets for flli 

' 2£ \l/2 



|t^) = -# T E/^(-) 1 

Pyx) ^— ' J m 



x 



{(x\Es(+))(Es(+Mr A {( x \Es(+))(Es(+)\m 
2tt 



P{x) 



J2 I dE^-fl\x\Es{+)){Es{+)\^ 



J-{ aTg ( x \E S (+))+ a Yg(E S (+)\iP)} (21) 
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III. THE ENTRANCE INTO THE INTERACTION REGION 



We start here to analyze the theoretical predictions of our formalism. To begin with, we consider the simple case 
of an initial Gaussian state prepared at t = in a zone where V(q) = 0, and directed towards the interaction region. 
This wave packet ip of width Aq — 26, is centered at qo < -well to the left of the onset of the interaction- with mean 
momentum po > 0. In configuration and momentum spaces we have: 

m = (^) 1/4 e- s2 Poe-^-^ (22) 

(pU) = (^)l/4 e -5 2 ( P -Po) 2 -* P9 o 
TT 

respectively. For appropriate values of qo,po and 5, such that pq5 » 1 and \qo\ >> 8, almost all the packet is initially 
at the left of the origin and moving with positive momentum towards the right. We use this simplifying assumption 
(the neglect of the Gaussian's tails with q > 0,p < 0) in our qualitative arguments, and in the intuitive descriptions 
of the processes that we will develope below. This will be indicated explicitly in the formulas by the use of ~ instead 
of =. However, we shall work with the full expressions ( p2[ ) wherever necessary in the calculations. For simplicity, we 
consider that the potential vanishes to the left of the origin. Preparing the state tp as s& id above with ip(q) « for 
q > 0, and its Fourier transform V>(p) « for p < 0, we have (Es(+)\ip) w S sr (j^) 1 ^ 4 ^p(p), so that 

(txa\rl>) m S sr J dE e- lEt (x\Er{+))ijj(p) (23) 

valid for the full range of values of x. Now, the initial state contributes to the time of arrival ( pi] ) a quantity 
d/dE a.Tg(Es(+)\ip) « —mqo/p, the same that in the free case. 



A. The quasi-classical case 



We start with the simple but illustrative case where the potential departs from for positive q with V(0) = 0, and is 
so smooth that the WKB method is valid. Then, for E > V(x) and to lowest order, one can neglect the exponentially 
small reflection that would vanish classically, getting energy eigenstates of the form 



(x\Er(+ 



»^<- i >v / ^ e "" +,,(I) v / ^j e ' /;j """ 



where p(q) = y/2m(E — V(q)). To this order and with a properly normalized wave packet as ours, (20) gives 

P 



P{x) w 6{-x) + 9{x)P + (x), P+(x) 



dp 



p(x) 



(24) 



(25) 



so that -f-r \ip(p)\ 2 is the (unnormalized) probability of arrival at the point x with momentum p{x). For the probability 

pyx ) 

in times of arrival one gets 

2 



P(t,x) 



6{-x) 



2tt 



6{x) 



2ttP + (x) 



dpe- M i/j(p) 
dE 



V p( x ) 



(26) 



which is the same as that of free particles for x < as corresponds to this order of approximation in which reflection 
is neglected, so that there is no information about V at the left of the origin. Finally, 



t a |V> / dp\^(p)\ 2 -{x-q } 



P+(x) Jo ^ P( x ) 



dp ■ 



H(p)\ 2 {- 



mq 
P 



x dq_ 

o p(qy 



(27) 



Therefore, for negative x we recover the toa of the free particle. What the above expression gives for x > is nothing 
else than the classical time of arrival at x, Eq. ([!]), for initial conditions (qo,p) weighted by the probability of these 
conditions. 
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B. Step potential and Hartman effect 



In general, the approximations that led to ( ]24| ) do not hold. For instance, reflection has to be taken into account, or 
V is such that the semiclassical approximation is no longer valid, etc. In any case, the particle may eventually reach a 
point q where E = V{q). Any further penetration beyond that point is a quantum fenomenon worth to investigate in 
terms of the toa. We address this question by considering a step potential V(q) = 8{q)V intercepting the path of the 
wave packet ip. We will then analyze the fate of the components of the wave packet with p > py = v2mV and with 
p < py. Classically, a particle in the first group will arrive with momentum p' = a/ \p 2 — P v | at the points x > 0, 
while one in the second group will bounce back at q — 0, without penetrating to the right. In the quantum case, one 
has for i<0a superposition of both, reflection and transmission, regardless of p/pv, while for x > one has 

(txa\ii>) » -fc /" dE{ — )V4 e - iBt x 

{9(E - V)T >e ip ' x + 8(V - E)T <e - p ' x }^(p) (28) 
where T> = 2p/(p + p') and T< = 2p/(p + ip'). Then, 

roc ppv 

P(»)«/ dp\T > iJj(p)\ 2 + dpe- 2 ^|T<^(p)| 2 (29) 

Jp v Jo 



is the probability of arrival at x, while 

P(t,x) 



2ttP(x) 



J dE(^)V*e- iE ^ 



{8{E - V)T >e ip ' x + 9{V - E)T < e- p ' x }i)(j>) 2 (30) 



jives the probability distribution in toa of the particles that arrive at this point. Finally 

1 



tx 



P(x 

rPV 



a \rr It M2r m Qo , mx , "irfarg(T>) 

dplTytpipJl { + — + } 

Pv P P P dp 

, -2 P 'x lrr 7/„m2; mq m darg(T<) 

dpe p T<V>(p) { 1 - A } 

p p dp 



(31) 



In the case of low potential steps pv « Po (c.f. Eq. (|22j)), where one can neglect the integrals over the interval 
[OjP^]) the probability of arrival reduces to the average of the transmision coefficient 1 | 2 , which is independent of 
x as corresponds to a transmitted free particle. T> is real in this case, so that ( - 0|t :c |'!/') is given by averaging over p 
the time spent to go from qo to at momentum p plus the time spent to go from to x at momentum p' . The only 
effect of the step is the reduction of the momentum from p to p' . 

In the opposite case where py >> po, only the integrals over [0,pv] give a sizeable contribution. The probability 
of arrival vanishes (exponentially) beyond the distance Ax = y associated through the uncertainty principle to the 
difference AE between the energy of the step and the energy of the particle. One then expects to detect a relative 
of this fenomenon in the time of arrival. In fact, the time spent from to x is given here through y d al J p T< = 
which is independent of the distance x, that is replaced by Ax. This is a case of the Hartman effect that here arises 
from the change 

P „iv' x P „ — v'x (32) 



p + p' p + ip' 

in the energy eigenstates as p crosses py from above. In short, the effect is a consequence of the fact that the phase 
is independent of x for p < py . 

In the general case one should take into account both contributions to (|3l]) . The relative importance of the second 
contribution in the rhs would depend on po — py and will always decrease exponentially with increasing x. However, 
a proper analysis of this situation calls for a description of particles better that that provided by first quantization 
and wave packets. We will defer this question to the next section where we discuss tunneling, the instance where the 
particle may reappear again beyond some point. 
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IV. ARRIVAL AT THE OTHER SIDE 



In this section we will study the modification of the times of arrival of quantum particles that traverse potential 
barriers. Our treatment deepens on the current understanding of the tunneling and dwell times. The literature is 
full of ad hoc heuristic arguments often disconnected from the standard mathematical and interpretative apparatus 
of quantum mechanics, whose value is therefore difficult to asses, as is their comparison with experiment. Here, we 
will follow the standard quantum mechanical treatment of Section 2. 

The time of arrival at a point x will now be given through a probability amplitude 



= / dE( — ) 1 / 4 e~ lEt (x\Es{+))(Es(+)\ij) 



We prepare the initial state as usual (as a right mover at the left of the barrier, c.f. above Eq. 
approximate (Es(+)\tp) w S rs (j^) 1 ^ 4: 'ip(p). The scattering state of relevance in ( |33| ) is given by 



(33) 

We again can 



(q\Er(+)) = 



9{-q){e*» + R{ P )e- ipq ) + 0(q)0(a - q)A(q,p) 



2np 

+ 6(q - a)T(p)e lpq 



(34) 



Expression valid for an arbitrary potential barrier contained in the range (0, a), where A(q,p) solves the appropriate 
Schrodinger equation with energy E — p 2 /2m. Also, T(p) and R(jp) are the transmission and reflection coefficients of 
the barrier. For a barrier of infinite range, the first and third terms in the rhs of (Q) should be better understood as 
asymptotic limits. 

Finally, in the case where x is at the right of the barrier, the amplitude can be approximately given by 



(txs\ip) 



/2vr 



dE(^)^e-« Et -^T(pMp) 



(35) 



The normalized probability density in times of arrival at x counts all the particles eventually recorded at x and only 
them, that is, the transmitted particles. According to it is given by 



(36) 



1 



2ttP(x) 

where we have normalized dividing by P(x), the total probability of arrival at x in whatever time t 



/+oo 
-OO 



dt |(tas|^)| 5 



dp\T(p)Tjj(p)\ 2 



(37) 



that is independent of x in cases like this, where x is beyond the range of the potential. In addition, it approximately 
simplifies to |T(po)| 2 for narrow wave packets with mean momentum po not too close (by above or by below) to the 
barrier momentum py = \/2mV . After a straightforward calculation we get for the average time of arrival at the 
other side of the barrier 



that can be written as 



2P(x) 



dE 



(—)V*j**T(p)fo) 



P&) Jo 



(38) 



d_ 

dE 



(_)l/4 e *x T(p )^(p) 



dp\T(p)i>{p)f^{x-q 



V 



dwg(T{p)) 
dp 



} 



(39) 



an expression thas has appeared before in the literature sometimes supported by heuristic arguments alone. It can 
be understood as the average value of the Wigner time [E2J over the transmitted state. 
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We will illustrate the predictions of the formalism for a simple square barrier of height V and width a. The 
transmission coefficient is in this case: 



T{p) 



2pp'e~ 



i pa 



2pp' cosp'a — i{p 2 + p' 2 ) sinp'a 



(40) 



where p 1 — \/p 2 — p v , that is imaginary for p below pv ■ Note the contribution —pa to the phase of T{p). This will 
substract a term a to the path length x — qo that appears in (|3S|). The barrier has effective zero width or, in other 
words, it is traversed instantaneously. This is the Hartman effect for barriers. To be precise, the effect is not complete, 
it is compensated by the other dependences in p'a present in the phase of T(p). In fact, it dissapears for (py /p) — > 0, 
where all the a dependences of the phase cancel out, as was to be expected because the barrier effectively vanishes in 
this limit. In the opposite case (p/pv) — > the effect saturates and there is an advance ^ in the time of arrival of 
transmitted plane waves, that turns into unexpected results for intermediate barrier momenta. 

We present our results for the time of arrival of the transmitted particles in Figs. 1 and 2. We consider the same 
initial state in both cases, namely the Gaussian wave packet of ( p2| ) with qo = — 30, po = 2,6 = 10, and m = 1, (we 
always use the natural units of the problem with Ti = 1). We have computed the time of arrival of the wave packet at 
x = 50 for an assortment of potential heights and widths, and have chosen the contents of those figures to highlight 
the most important results. 



C5 
# > 

'u 

•— 

a 
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O 
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FIG. 1. Average time of arrival at the other side of a barrier of fixed width a — 15 as a function of the barrier momentum 
pv- The parameters of the initial Gaussian wave packet are go = — 30, po = 2, S = 10, m = 1, and the arrival at x — 50, in 
units with h — 1. The solid line is the quantum average (pa), the dashed line is the phase time with momentum po. The 
asymptote to the left (pv — > 0) is the time of arrival for free particles (to = 40), the one to the right is the Hartman time 
t H (a) = to - = 32.5. 



We show the time of arrival at the other side of a barrier of momentum py in the range a — (1.6, 2.6) in Fig. 1. 
For incident plane waves with momentum po, the barrier would be crossed over for pv < po, and tunneled through 
for py > pq. Some retardation would be expected in the first case, just because the travel over the barrier would be 
slowlier than the free travel. This is clearly seen at the left of po in the figure. Classically, the delay would grow from 
zero (time to) to infinity as py grows from to pq. The quantum behaviour is similar, with the oscillations of the 
phase time swept away by the average that remains finite. To the right of po, there is a dramatic difference between 
the Wigner result, that inmediately sticks to the Hartman prediction tu, and the wave packet result, for which the 
time continues to increase up to a certain barrier height and then, suddenly, drops to tn ■ This strange behaviour can 







be explained in the following manner: Not being monoenergetic, the wave packet has momentum components above 
and below py. The first of these cross above the barrier, get retarded, and are responsible for the high time value for 
Pv just to the right of po- However, as the barrier continues to grow, they become an ever lesser part of the packet. 
The other parts of the packet (the components with momentum p < py) tunnel through the barrier, and experience 
the Hartman advance. They would arrive at x in a time tn- Their relative importance in the wave packet increases 
steadily as py continues to grow and, eventually, they overcome the retarded components and the process becomes 
pure tunneling. Then, the time of arrival drops to tn- We have numerically checked this behaviour, that we have 
analyzed for several values of the barrier width in the range (2,30). All the results are similar: Monotonic grow of the 
time from py = (where t — t ), up to py w 2.5, where t drops suddenly to tn- The general trend is a slow increase 
in the value of the barrier momentum py at which the drop takes place, that shifts from about 2.2 to 2.7 as a changes 
from 10 to 30. The maximum value of the time of arrival (t x ) that is obtained just before the drop also increases; it 
is around 95 for a = 10 and around 450 for a = 20. 




2 4 6 8 10 12 14 

Barrier thickness 



FIG. 2. Average times of arrival at the other side of a barrier as a function of the barrier width a. The initial wave packet 
is the same as in Fig.!, and again x = 50. We show the predictions for pv = 0, 1.6, 1.8, 2.2, 2.3, 2.4, 2.6. The solid lines are the 
quantum averages (p9f), the dashed lines are the corresponding Wigner times for a particle with momentum po 

We show in Fig. 2 the average time of arrival and the Wigner (phase) time as a function of the barrier width a in the 
range a — (0, 15). We display the predictions for different barrier heights py = 0, 1.6, 1.8, 2.2, 2.3, 2.4 and 2.6. For the 
free case (py = 0, or a = 0) all the results converge to to = 40. We now disscuss the solid lines (t x ). The oscillatory 
curves above to correspond to py < po- The get steeper as their momenta approach po from below. The curves that 
stand partially below to correspond to py > po (tunneling). They share a similar behaviour: As the barrier width 
grows from a — the time of arrival decreases, practically saturating the Hartman time t(a) ~ tn(a) = to — 
Then suddenly, at a certain width (that increases with py), the average time jumps dramatically to values that 
correspond to a long retardation. Note that the jump for py — 2.6 lies outside the range of the figure. This 
behaviour is complementary to that shown in Fig. 1. Here, for py > po and moderate a, tunneling is the dominant 
phenomenon and the time average tends to reproduce tn- However, as the barrier gets wider, tunneling gets more 
and more depressed. In comparison, the intensity of the retarded components that pass over the barrier is basically 
independent of a. They get relatively more and more important and, eventually overcome tunneling, giving rise to the 
observed transition. In practice, for wide enough barriers, the probability of tunneling vanishes, and the other side 
can be reached only by the very improbable and very slow travel over the barrier. This behaviour has been noticed 
independently in ]33| ], and explained in the same way. In addition, we have the tools to check these explanations. 
In particular, the first product of our formalism is P(t,x), the probability distribution in times of arrival at x. Our 
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numerical analysis for x = 50 and the different py's and a's that we are discussing here show similar almost Gaussian 
shapes for these distributions, as correspond to the initial wave packets chosen, and similar widths for these P(t,x), 
whose maxima are placed close to the corresponding mean values (t x ). As expected, the probabilities get numerically 
smaller as the corresponding events become more and more unlikely. In short, these distributions give the best support 
for the validity of the explanation offered here for this striking behaviour, that can be understood only after weighting 
the obtained time of arrival with the relative probability of the actual event to which it corresponds. 



V. QUANTUM REFLECTIONS 



Having analyzed the modifications introduced by the transmission phenomena in the time of arrival at the other 
side of potential barriers, we turn to the case of reflection. We divide the analysis into the two seemingly different 
cases in which there is classical reflection, and in which it is absent. The first case is characterized by the presence 
of at least one turning point in the path of the particle. The second one, by the absence of any of them. Quantum 
mechanically there could be some transmission in the first case, and some reflection in the second one. Accordingly, 
we separate the disscussion that follows into the two main disjoint cases that cover all the possibilities. These are 
the case where the potential energy grows to infinity somewhere (total reflection), and the case where it is bounded 
everywhere (with partial reflection and transmission). 



A. The case of total reflection 



The potential energy could grow unbound, thus reflecting any conceivable incoming state. We consider here a 
monotonic potential energy that vanishes for q — > — oo and goes to infinity for q — > oo so that \im q —, +00 (q\E) = 0. 
This removes the degeneracy of the energy eigenstates. As no state may arrive from the right, (q\El(+)) = 0. 
The eigenstates \E) will contain the same amount of positive and negative momenta, so that their asymptotic form 
normalized to one traveling particle per unit time is lim g ^_ 00 (<7|£') = -^(^g.) 1 / 4 cos(pq + 5(E)), where 5(E) is the 

phase shift. This also fixes completely the eigenstates for finite values of q. 
The time of arrival at an arbitrary point x is now 

(1>\t*\fl>) = ^JdEyJ&\(x\E) W)| 2 X 

^{arg^+arg^}} (41) 

which is the average of a quantity independent of x\ This comes about because in the present situation the reflection 
coefficient R = exp(— 2i5) is unimodular. Then, the net current density vanishes, so that &rg(x\E) is independent 
of x. This is the quantum version of the classical result that the sum of the times of arrival at x of the incoming 
and returning particles is twice the toa at the turning point and so, independent of x. Obviously this ceases when 
\R\ becomes smaller than 1 (so that the net current density is finite), something that is possible only when V is 
finite everywhere. Even then, the classical result is recovered from the quantum case in the limit (E/V) << 1 where 

\R\ - 1. " 

The individual times of arrival of the incoming and the returning particles can be obtained straightforwardly by 
writing the enegy eigenstates as 

(q\E) = 4^) 1/4 M(q,E) cos<f>(q,E) (42) 



where M is a real function with lim g ^_oo M(q, E) = 1, that vanishes faster than an exponential for q — > +oo to satisfy 
the asymptotic form of the Schrodinger equation. The state is thus written as the superposition at each point of an 
incoming and a reflected wave with equal amplitudes, so that the net current vanishes everywhere. The phase cf> is 
fixed by lim g ^_ 00 <fi(q, E) = pq + 5(E) to match the asymptotic form of the eigenstate disscussed above. Its derivative 

gives the two opposite velocity fields v± (q, E) — ± d< ^^ interfering at q. We recall that this exact expression is valid 
for all the potentials of the form we are considering here. The probability of ever arriving at x and the toa can be 
given by straightforward application of ( pp[ ) and ( |2l| ) by 

P(x) = J dE M 2 (x,E) cos 2 <f)(x,E)\{E\ip)\ 2 (43) 
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= 2p^y J dEM 2 (x,E) cos 2 <t>(x,E) x 

\(e\i>) | 2 [tifo £)+*»•(*>£)] 

which is the weighted average over energies of the times of arrival of the incoming and the reflected waves: 

ti(x, E) = ^{<f>(x, E) + arg(£|V>} 
tr(x, E) = -^{~4>{x, E) + arg(£?|V)} 



(44) 

(45) 
(46) 



whose sum is explicitly x independent. 

To illustrate these results we consider now the case of a potential that vanishes at the left of the origin and is linear 
at the right, i.e. V(q) = 0(q)fq, where / is the force exerted on the particle. This could be a model for a (charged) 
particle in a constant electric field, or in the gravity field of the Earth. In this case one gets M and <j> m terms of the 
Airy function Ai and its derivativeAi'. 



M(q,E) 




Ai[ z ]2+(-f)2Ai'M 



1 for q < 

2 

for q > 



(47) 



Ai[z ]2+(^)2Ai'[ 20 ] 2 

where z = k/q — p 2 /k 2 ,zo = —p 2 /k 2 with kf — (2m/) 1 / 3 . For the phase one has 



,E) 



arctan^^^TR ) F,,r " 



A ' 

arctan(— fc/ ^jTj ) for g > 
so the phase shift is given simply by 6(E) — 0(0, E). 



(48) 



a 



ca 
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u 




. 5 

Time of arrival 

FIG. 3. Probability distribution P(t,x) in times of arrival at x, for x = go (solid line), x — 0.5 go (dashed line) and at the 
classical turning point x — E/f (dot-dashed line). The vertical lines correspond to t = and to the classical turning time 
respectively. The distributions are bimodal, with the two peaks corresponding to incidence and reflection getting closer as x 
approaches the turning point. 
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We present in Figures 3 and 4 our results for the the case of a force of nominal value / = 100, being the parameters of 
the initial Gaussian wave packet ( |22| ) qo = — 2,po = 10, 5=1 and m — 1. For the normalized probability distributions 
in times of arrival ( |l9| ) we get pairs of peaks of equal heights - as correspond to total reflection - that tend to merge 
into one as the detector is displaced towards the classical turning point. This behaviour of the peaks is also observed 
for the averaged times of arrival, that follow the classical times. The small deviations from the parabolic form are 
negligible in comparison with the widths of the distributions shown in Fig. 3. We have explored numerically the 
details that change uninterestingly according to the values of f,po,S etc. so, we do not show them here. The general 
picture is always the same: at the far left (|go| >> E/f) the potential acts as an infinite height wall. The only 
sizeable consequences of the actual strenght of the force are felt at positions between the origin and the turning point, 
where they resemble the classical effects. Part of this comes from the fact that here position and energy combine into 
only a variable q — Ej f. But the resemblance arises because total reflection is always present here, quantum as well 
as classically. This will be more clear in the next section where we consider partial reflection that lacks of classical 
analog. 




-2 -1.5 -1 -0.5 0.5 

Detector position 



FIG. 4. Average time of arrival for different detector positions. The vertical lines correspond to x — qo, x — and x — E/f 
respectively. The solid line is the classical time, the dashed line is the quantum average of i* (lower part) and of t r (upper 
part). 



B. Partial reflections 

In classical mechanics a potential interaction energy speeds up or slows down the particles according to the local 
value of the force F(q) — — dV Q^ ■ Accelerated or decelerated, the particles continue to move along the same path 
without reversing the direction. Only when one of them intercepts a turning point (i.e. a point q where E = V(q)) 
the particle bounces back or, in other words, is reflected with probability Pr = 1. In the absence of these points, the 
particle is always transmitted with probability Pt = 1. Thus, most of the time Pr = 1, Pr = 0. Only at the turning 
points Pt = 0, Pr = 1. 

Quantum dynamics offers a very different perspective of the motion of the particles. The Schrodinger equation 
implies that at every point where the potential energy is finite, the particle is partially transmitted and partially 
reflected, that is < Pr < 1,0 < Pr < 1, with Pt + Pr = 1. The case of total reflection analyzed in the previous 
section is one of close correspondence between the classical and the quantum results, as we shown there. Interesting 
departures from the classical behaviour arise when there is no classical reflection. We will analyze this case here. 
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To fix ideas, we consider a well behaved potential energy V(q) > V finite q, that vanishes at the spatial infinity 
faster than q^ 1 . In these conditions the energy eigenstates can be written everywhere as a well defined superposition of 
transmitted $t r (<7, E) and reflected $ re /(g, E) waves, characterized by the positive or negative value of their currents: 

~2^( $ tr ^ ®tr) > 0, and -5S:($* e/ 25 $re/) < 0, with different amplitudes |$ tr | ^ |$ re /| as corresponds to this 
case of partial reflection. The eigenstates of interest can be written as 



/ TYl 

(q\Er)=^—{$ tr (q,E) + <S> ref (q,E)} (49) 

These waves are univocally determined by their asymptotic conditions, namely: 

lim $ tr (q,E)= e lpq , lim $ tr (q, E) = T{E)e lpq 

q — > — oo q— ' + 00 

lim <P ref {q,E) = R(E)e-^ pq+25 ^\ lim $ ref (q, E) = (50) 

q — > — oo q — >+oo 



as is the case for an incoming rightmover (49). The results of the previous section arc recovered in the limit where 
T(E) — ► which is the case only if the potential energy grows to infinity somewhere. 

If we prepare our initial Gaussian state ip(q) at a point q — qo where the potential energy is smooth enough, and 
keep the initial momentum po > large enough to consider ip(p) ~ for p < 0, we can use the approximations 

(Esty}n6 rs J^$t(qo,E) $(p)| «* r .Jje-*« \Mp)\ (51) 

We have used the second of these already in Eq. (g^) . It is valid when V(q) ~ for q in the qo neighbourhood where 
ip(q) is sizeable. We assume this is the case in what follows. 

One of the deepest consequences of the superposition of transmitted and reflected components that makes up the 
eigenstate ( |49"| ) is that it leads to the inescapable presence of interferences. In fact, the probability of presence at a 
point q, and other quantities depending on it, contain the sum |(g|i?r)| 2 cx \&tr\ 2 + \®ref\ 2 + 23?($ tr <I>* e j), whose 
last term is the interference term. One could say that, everywhere in its motion through the interaction region, the 
quantum particle will be found in an evolving entangled state of transmitted and reflected components. This can be 
traced back mathematically to the continuity of the solutions of the Schrodinger equation and of their first derivatives, 
and to the associated Wronskian theorem. Physically, this may introduce all sorts of interpretative difficulties in the 
analysis of particle motion. 

Summarizing, interferences pervade the realm of quantum motion. They will show up in almost every quantum 
mechanical situation. Our analysis of the time of arrival is not an exception. We have avoided refering to them till 
now by focusing on very specific cases. These were: The choice in Sect. 3.1. of a very smooth potential analyzable 
semiclassically by the WKB method, that neglects reflection. The analysis in Sect. 3.3. of the time of arrival at 
points located at the other side of the barrier, where $ r e/ = so that any interference with the transmitted wave 
vanishes. Finally, the analysis made in the previous section, where we just ignored the effects due to the overlap of 
incoming and reflected waves in P(t,x), and the lack of a clear cut separation between t% and t r in the presence of 
interferences. To be precise, we dealt with reflection without paying the due attention to these subtleties. We repair 
the ommission here. 

The amplitude in time of arrival at a position x within the interaction range can be given by using (B^) and j5l|) 



in (33) 



(txs\ijj) = {A tr (t, x) + A ref (t, x)} 



'2tt 



dp J |^(p)| e -'(«*+P*) { $ tr(Xj E) + $ re/(Xj E)} (52 ) 



This gives for the probability of ever arriving at x Eq. (|20|) the sum of three terms: The two separated probabilities 
Ptr, P r ef of arriving with positive or with negative current density, and a quantum interference term, whose presence 
deprives the previous two of direct physical meaning. We thus get P(x) — Pt r (x) + P re f(x) + I{x) with 



P tr (x) = dt\A tr (t,x)\ 2 w / dp m P )\ 2 \<S> tr (x,E)\ 2 (53) 

re/ J re/ J re/ 



and an interference term 
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I(x) = 2 J dt$t{A tr (t,x)A* ef (t,x)} 

« 2 J dp \^(p)\ 2 ^{e- ipc > $ tr (x,E) <S> ref (x,E)*} (54) 

The above quantities depend on the probabilities of transmission or reflection from the initial position qo to the 
actual value x. Consider a bounded potential barrier of finite range, but otherwise arbitrary. Behind the barrier P re f 
vanishes, while P tr is given by fl37| ) with a value independent of x, but strongly dependent of p ,S and of the barrier's 
height and width. For x at the left of the barrier $ tr = e lpx (what we are denoting as transmission is here incidence), 
but $ re/ = R(E) e- l{px+25 ^\ and only when there is no reflection (no barrier) the intereferences dissappear. For 
the total reflection case of the previous section, we get Pt r — Pref, while the interference term gives rise to the term 
cos{2(px + 5(E))} that builds up the factor cos 2 <fi that appears in (43) and (|44|). However it does not prevent the 
definition of the quantities ( |45| ) and(f46[) that allowed to split the toa (44) into two positive contributions interpretable 
as the independent (t x ) of an incoming packet and a reflected one (Fig. 4). 
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FIG. 5. Probability distribution P(t,x) in times of arrival for reflection from finite potential barriers. The initial wave packet 
is the Gaussian one with m = l,po = 2, 5 = 10, placed at go = —150. The arrival position is at x — —100. The upper left 
figure is for a barrier width a — 4 and pv = 2.2. At its right is the case a = 4 and pv = 1.9, that is enlarged in the lower 
left part for the range t = [100, 180]. An illustrative case of multimodal reflection distribution is shown at the lower right part, 
that corresponds to pv = 1.9 and a = 6. The vertical grid lines correspond to the phase times of the incident wave e lpax and 
of the reflected waves: e~ rp ° x for pv > po, and the superposition sinp a e~ tp " x for pv < po- 



For finite barriers reflection is always present with an energy dependent coefficient R(E) < 1; it is less probable 
than incidence, and tends to vanish as the barrier does. In Fig. 5 we give the probability distributions of toa P(t, x) 
at a point x, whose bumps indicate, as in Fig. 3, the arrival of incident and reflected parts of the time evolved initial 
wave packet. This is the Gaussian one with m — l,po — 2.5 — 10, placed at qo = —150. The arrival position is at 
x = —100, far from q$ to avoid interferences. The two upper figures are for a barrier of width a = 4. At the left is 
the case where pv — 2.2, and at the right that with pv = 1.9. In both cases there is an incidence bump centered 
at t = m(x — qa)/po = 25, and a structure to its right corresponding to reflection. For py — 2.2, and for all the 
cases of total classical reflection (pv > Po), the latter is a Gaussian- like bump shifted from the classical value at 
t = m{—x — qo)/po = 125 by an amount (y(^))- However, for py = 1.9 (in general for py < Po), the reflected 
distribution has a multi-bump shape difficult to understand in terms of the phase time or of any other approximation. 
In particular, neither the number of peaks, nor their positions heights and widths can be approximated by straight 
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stationary phase methods. Two illustrative cases of these shapes are shown in some detail in the two examples of the 
lower part that correspond to py = 1.9 and two close widths a = 4 and a = 6. 



VI. CONCLUSIONS 



We have worked out a formalism for obtaining the time of arrival at a space point of particles that move through 
interacting media. Our construction follows a circuitous path: we desist from first computing the classical toa of the 
problem, and then quantizing it, a procedure that leads to a dead end. Instead, we start from the quantum toa of 
the free moving particle, and then transform it canonically to the interacting case. This is achieved by the use of the 
appropriate Moller operator that implements the quantum version of the Jacobi-Lie canonical transformation to free 
translation coordinates in phase space. In the classical case we have the transformation of Eq. (|5^) whose quantum 
counterpart is 



where Ho and H. are the Hilbert spaces of the free and interacting particles, and Hq,H the respective Hamiltonian 
operators in these spaces. For simplicity, we have only addressed explicitly cases in which the transformations are 
unitary, which is the case when cr(H) = ct(Hq). More general situations that require of isometric transformations, 
deserve a separate treatment by their physical relevance. 

What we obtained here is a quantum formalism for the toa in terms of a POVM given by 



which measures the probability of arrival at x during the time interval (ti, £2). The normalized probability distribution 
P(t,x) was given in the Eq. ( |T9| ) of Sect. II. B. Our results are thus within the standard formalism of quantum 
mechanics and can be interpreted in the standard way. There is nothing special that singles out our theoretical 
predictions as unsuitable for comparison with the experimental results. On the contrary, our formalism predicts the 
result of actual experiments in the form of numeric values and statistics for the recorded events. 

After the definition and theoretical analysis of Sect. II. we have performed explicit and complete calculations for 
the cases of an unbounded linear potential, of the step potential and of the square barrier. Our analysis of the 
quasi-classical case shows that in this limit the toa is simply given by the average of the classical time of Eq. (jl]) 
over the quasi-classical wave function. In the case of reflection, and for the arrival point placed between the initial 
position of the wave packet and the turning point (x < 0), the probability distribution P(t,x) is governed by the 
quantum superposition of the incident (A tr ) and the reflected (A re f) wave packets. In the case of total reflection, 
where both are equally probable Pt r (x) = P re f(x), we have obtained separate positive (t x ) even when both amplitudes 
overlap. These were interpreted as the toa's of the incident and reflected particles, and compared successfully with 
the classical prediction. For partial reflection, P re f(x) < Pt r (x) non overlapping amplitudes are necessary to to get 
separate average values for these times. This problem is shared with the position and other operators. It is not a defect 
of the formalism, but an effect of the interferences. Fortunately enough, our formalism provides us with the probability 
distribution P(t,x) whose diverse humpy-bumpy shapes (Figs. 3 and 5) give the most complete information of the 
posible experimental outcomes. 

In the course of our numerical analysis we have detected that the phase time not always gives a good approx- 
imation to the most probable time of arrival. It provides a first estimate of the time spent in the transmission or 
reflection, after substracting the time of free flight. For transmitted wave packets we have reobtained the advancement 
(i.e. a decrease in the toa) in the case of pure quantum tunneling. This phenomenon, predicted by Hartman long 
time ago I^H , has been experimentally evinced by the two photon experiments at Berkeley [|l8],[l9| and the tunneling 
of optical pulses at Wicn pQ] , However, our formalism predicts a striking departure from the Hartman bound that 
we explain in detail in Sect. IV. Our results for square barriers neatly show the expected advancement roughly pro- 
portional to the width At — —ma/p (Figs. 1 and 2). However, whatever the mean energy (E < V) of the incident 
wave packet, there is always a width ao such that for a > ao the (very retarded) components of the packet that 
stand above the barrier dominate over the (probabilistically very depressed) tunneled ones, giving an overall effective 
strong retardation. In other words, when the barrier is wide enough, its width dominates over the Hartman lenght Ax 
disscussed above Eq. ([32]), that has a purely quantum origin. This restores the classical expectation of no tunneling 
and very long delays. 

We have also found other unanticipated phenomenon for purely quantum reflection: the multiple bump structure 
that appears when py < pq. We have shown in Fig. 5 this structure, that in some sense is a counterpart of the 



{H,H} 




{n ,H } 



(55) 




(56) 
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interference pattern that appears in multiple reflection of stationary waves. We think that this feature, even if 
less spectacular than the superluminal tunneling of photons, deserves experimental confirmation. An appropriate 
modification of the two photon experiments could serve for this purpose. It would require to place a quantum mirror 
in the path of one of the entangled photons, and check for the presence (or absence) of the multiple dip structure in 
the number of coincidence counts predicted by the formalism. 

All the examples above show that our construction of a quantum toa operator suitable for the presence of interactions 
allows the exploration of many physical details in relevant situations. Its extension to higher dimensional cases poses 
no conceptual difficulties and opens the possibility of treating new questions. Of great theoretical and experimental 
interest will be the extension of this formalism to the cases in which the Hamiltonian has bound states, where isometric 
(instead of simply unitary) transformations will be requiered. 



APPENDIX A 

In the modern literature J37J , a classical Hamiltonian system with n degrees of freedom is called completely integrable 
( a la liouville) when it satisfies the conditions a and b below: 

a. There are n compatible conservation laws $>i(qi, . . . ,q n ,pi . . . ,p n ;t) = Ci, i = 1, . . . ,n, that is: 

a.l. $ i = {$ i ,ff}+^i = 0, Vt = l,...,n. 

a. 2. {$ i ,$ i } = 0, Vi,i = l,...,n. 

b. The conservation laws define n isolating integrals that can be written as: 

b. l. <J>i = Ci => Pi = <f>i(qi, . . . ,g„,Ci, . . . , C n ;t), V i = 1, . . . ,n. 

d4>j _ d<pj 
dqj dqi 



b2 Vi i-l n 



In these conditions, Hamilton equations define an integrable flow, that is, a system of holonomic coordinates (q(t),p(t)) 
in phase space for each instant of time: 

Qi(t) = qi(qo,Po;t), i = l,...,n. (57) 
Pi{t) =Pi(qo,Po;t), i = l,...,n. 

In other words, given a set of initial conditions (qo,po) of the system, at each instant of time t the system arrives at a 
point {q{t),p{t)) in phase space. Conversely, these points define the corresponding times of arrival. In this case, time 
meets the requirements to qualify as a derived variable in phase space. 

As Lie pointed out, for any arbitrary time there is a special choice of coordinates in phase space that mathematically 
eliminates the effects of interactions from these integrable flows, (the new positions are ignorable coordinates). More 
simply, integrable systems are canonically equivalent to a set of translations (or circular motions) at constant speed. 
It is customary to denote the variables that determine these translations as action- angle variables, which strictly is 
appropriate only in the case of periodic systems, where the (closed) flow lines are topologically equivalent to circles. 

For integrable flows, there is a canonical transformation (the Jacobi-Lie transformation) 

{q,p; H(q,p)} {Q, P; H(Q, P)} (58) 

with H(q,p) = H(Q, P), that gives the free translation coordinates P(t) = P, and Q(t) = —t+Q of the translation flow 

with H{Q, P) = j—, in terms of the coordinates and momenta (q(t),p(t)) of the actual flow with H(q,p) — 2^ + V(q)- 
This transformation is of the form W(q,P), that is, a function of the old coordinates and the new momenta, so that 

n dW dW 

Q= W P= ^q- (59) 

Finally, W can be obtained explicitly as a complete integral of the Hamilton- Jacobi equation: 

dW P 2 
oq 2m 

Now, the canonical relation among the new and the old variables is: 
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P = sign(p) y/2m H(q,p) 

Q = « 



1 _ Yla. 



where Qo is a constant. As a byproduct, time gets defined in equivalent manner in terms of the old variables, or of 
the new ones. If the particle arrives at q(t) = x in the instant t{x) = t, then: 

t(x) = ^(X-Q) = sign(p) f mdq ' (61) 

where X — dW(x, P)/dP (obviously, X — Q(t) by construction). This duality, devoid of practical interest in the 
classical domain, is at the foundations of the quantum method developed in this paper. Finally, note that for simplicity 
we have specialized the notation to the case of autonomous Hamiltonian systems with only one degree of freedom, all 
of them trivially integrable (H(q,p) = E being the needed conserved quantity). 



APPENDIX B 

For free particles Eq. ([!]) gives t x o(q,p) = m(x — q)/p that, in spite of its simplicity, presents some problems for 
quantization |3] ]l4]Jl6|| whose solution we outline here. First of all, it requires symmetrization: 

Mq.p) = m(- - i{q, -}+) = - e -P- & (62) 
P 2 p V P V P 

As is well known, the eigenstates |ixsO) of this operator in the momentum representation can be given as (h = 1) 

( P \txsO) = 6{sp)\ M exp(i£-t) (p\x) (63) 
V m 2m 

where we use s — r for right-movers (p > 0), and s — I for left movers (p < 0.) The label stands for free case. 
Finally, the argument sp of the step function that appears in the momentum representation is +p for s = r, and — p 
for s = I. The degeneracy of the energy with respect to the sign of the moment is explicitly shown by means of the 
label s in the energy representation, where 

(Es'0\txs0) = S s , s (— y/ 4 e lEt (EsO\x) (64) 
m 

Summarizing, there is a time (of arrival at x) representation spanned by the eigenstates 

\txaO) = (^)VV"°'n s0 |x) (65) 
m 

where II s o projects on the subspace of right-movers (s — r), or left-movers (s — I), i.e. 

/>oo 

n s0 = / dE\Es0)(Es0\ (66) 



These time eigenstates are not orthogonal, which in the past gave rise to serious doubts about their physical meaning. 
The origin of this problem can be traced back to the fact that ( |62] ) is not self-adjoint, that is (tp\t X Qip) ^ (t x0 (p\ip). 
This was proved by Pauli (T) long time ago and is due to the lower bound on the ener gy spectrum. The problem 
emerges as soon as one attempts integration by parts in the energy representation. Ref. |27| is a recent illuminating 
review of these and other related questions. 

The measurement problem posed by this not self-adjoint toa operator can be solved by interpreting it in terms of 
a Positive Operator Valued Measure (POVM), that only requires the hermiticity of t^o (i-e. t^o = (txo)* T )- Here, 
instead of a Projector Valued spectral decomposition of the identity operator, one has the POVM 

P (n.(xy,t u t 2 ) = dt \txs0) (txs0\ (67) 



V / 2 ^(^)V4 e ^o tns0 n(*)iTo 



e -iH t ^ 2H Q y/4 
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where n(x) = |a;) (x\ is the projector on x. Here, Po(l,2) 2 ^ Pq(1,2) because \txsO) (txsO\ is not a projector, as the 
states are not orthogonal, but where the limit as t — > oo of Pq(— t, +t) is the identity. The attained time operator is no 
longer sh arp, b ut is well suited for measurement. This solution has been implemented in and extensively analyzed 
in refs. J|4|,|35|] and in the review ]27|] . In this POVM formulation the toa is given by the spectral decomposition 

t (tfo,n(z))= f + °° dtt(^) 1 / 4 e lHat V (x)e-^ t ('^) 1 / 4 (68) 

■/-oo m m 

where Vo(x) = XL ^so n(x) n s o, which is not a projector. 
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